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Torsional  Oscillations  of  a  Gun  Barrel 


ABSTRACT: 

This  report  is  concerned  with  an  analysis  of  the  torsional  oscillations 
of  an  artillery  barrel  when  subjected  to  the  moment  of  an  accelerating 
projectile. 

A  modal  series  is  obtained,  the  coefficients  of  which  are  the  time 
dependent  generalized  coordinates.  For  purposes  of  obtaining  numerical 
results,  this  series  was  truncated  at  two  terms;  the  subsequent  system  of 
ordinary  differential  equations  was  programmed  for  the  analog  computer. 
Graphic  solutions  are  given  for  several  values  of  the  barrel  size  to  projectile 


mass  ratio. 
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NOTATION 

Coordinate  denoting  position  along  the  barrel. 

Coordinate  denoting  radial  distance  from  the  center  of  the 
barrel  cross  section. 

Angular  displacement  of  barrel  cross  section  and  coordinate 
direction  in  a  cross  section  perpendicular  to  the  radial 
direction. 

Time. 

Computer  time  variable  =  4000  t. 

Distance  of  the  projectile  along  the  barrel. 

Total  angle  of  rotation  of  the  projectile. 

Length  of  the  barrel. 

Shearing  modulus  for  barrel  material. 

Mass  density  of  barrel  material. 

G/p 

Polar  moment  of  inertia  of  barrel  cross  section. 

Mass  moment  of  inertia  of  the  projectile. 

Constant  determined  by  the  helix  angle  of  the  barrel  rifling. 
Torsional  spring  constants. 

An  appropriate  average  angular  acceleration  of  an  infinitesimal 
length  of  barrel  element. 

Shearing  strain  between  the  x  and  4>  directions. 

Shearing  stress  on  a  cross  section. 

Time  dependent  moment  exerted  on  the  barrel. 

The  n*"^  eigenvalue, 
th 

The  n  natural  frequency. 


normal  mode  of  vibration,  or  n*'^ 

Qjj  (t)  The 


eigenfunction. 


generlalized  coordinate. 


Chapter  I 


THE  EQUATION  OF  MOTION 


1.  1  Introduction 

The  analysis  of  a  complex  physical  system  inevitably  necessitates  the  use 
of  a  simplified  mathematical  model. 

The  physical  dynamics  problem  to  which  this  investigation  has  been  addressed 
is  the  problem  of  determining  the  torsional  oscillatory  behavior  of  a  rifled  artillery 
barrel  when  subjected  to  the  time  dependent  moving  torque  of  an  accelerating 
projectile. 

The  analogous  mathematical  problem  treated  here  is  the  problem  of  determin¬ 
ing  the  angular  displacements  of  a  cylindrical  shaft  when  subjected  to  a  time  depen¬ 
dent  line  moment  which  is  dependent  in  part  upon  the  angular  displacements.  The 
shaft  is  assumed  to  be.  made  of  a  homogeneous  isotropic  material,  and  the  deforma¬ 
tions  are  required  to  be  small  so  that  the  material  behaves  in  a  linearly  elastic 
manner. 

The  eigenfunction  technique  is  employed  in  that  solutions  are  sought  each 
of  which  can  be  represented  by  the  product  of  a  normal  mode  of  vibration,  a 
function  of  distance  alone,  and  a  generalized  coordinate,  a  function  of  time  alone. 

Two  systems  are  considered:  that  for  which  the  shaft  is  fixed  at  one  end  and 
otherwise  free;  and  that  for  which  the  shaft  is  free  at  both  ends,  but  elastically 
supported  at  two  intermediate  points  along  its  length. 

The  solution  of  the  problem  is  expressed  in  the  form  of  an  integral  equation. 
For  the  case  of  the  shaft  fixed  at  one  end,  the  modal  series  resulting  from  the  use 
of  the  eigenfunction  technique  is  truncated  at  two  terms,  and  the  resultant  pair  of 
second  order,  ordinary  differential  equations  with  variable  coefficients  is  solved 
through  the  use  of  an  analog  computer.  Numerical  results  corresponding  to  some 
typical  artillery  data  are  obtained  and  presented  graphically. 


2. 


1.  2  Derivation  of  the  Equation  of  Motion 

A  differential  equation  of  motion  for  a  circular  shaft  in  torsional 
oscillation  about  its  longitudinal  centroidal  axis  will  be  derived.  The  shaft  is 
assumed  to  be  made  of  an  elastic,  homogeneous,  isotropic  material.  The  angle 
of  twist  must  be  small  in  order  that  the  proportional  limit  of  the  material  not 
be  exceeded.  The  shearing  stress  on  any  transverse  cross  section  is  assumed 
to  vary  directly  as  the  distance  from  the  center  of  the  section;  this  requires  that 
plane  cross  sections  remain  plane  and  that  any  straight  line  segment  in  a  cross 
section  before  deformation  remain  a  straight  line  segment  after  deformation. 


Figure  1  shows  a  portion  of  this  shaft.  L  (x,  t )  is  the  torque  per  unit 
length  applied  to  the  shaft,  the  positive  direction  taken  counterclockwise  as 
shown.  ^  (x,  t  )  is  die  angle  through  which  a  particular  cross  section  rotates 
as  a  rigid  unit.  A  segment  of  the  shaft  A  x  units  long  is  shown  in  Figure  2. 

T  (x,  t )  is  the  internal  moment  exerted  on  a  cross  section  by  the  adjoining 


section  . 
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L (x,  t ) 


X,  t ) 


Figure  2 


The  sum  of  all  moments  on  the  element  of  Figure  2  must  be  equal  to 

*_  ♦ 

1  a  ,  where  I  is  the  mass  moment  of  inertia  of  the  element  about  its  centroidal 
axis  and  a  is  an  appropriate  average  angular  acceleration.  Thus, 

L  (x,  t  )  Ax  +T  (x  +  Ax,  t  )  -  T  (x,  t  )  =  I*  .  (1) 

d  t 

Expanding  T  (x  +  Ax,  t  )  in  a  Taylor  Series  about  the  point  (x,  t ),  we 
obtain 

T  (x  +  Ax,  t  )  =  T  (x,  t  )  +  (x,  t )  Ax  +  .  .  .  O  Ax  .  (2) 


Substituting  (2)  into  (1)  and  passing  to  the  limit  results  in 


I  /  \  9  T  ,  , 

L  (x,  t  )  +  (x,  t  ) 


pj  ^  (X,  t), 

9  t 


(3) 


where  p  is  the  mass  density  of  the  material  and  J  is  the  polar  moment  of 
inertia  of  the  cross  section  about  the  centroidal  axis. 

The  shearing  strain,  ,  is  defined  to  be  the  angle  (  *  Q  ),  where 

n  is  the  angFb  in  the  deformed  state  between  a  line  originally  parallel  to  the 
x-axis  and  a  line  in  the  <t>-direction. 


4. 


F igure  3 


For  small  deformations  (  ^  )  is  approximately  equal  to  tan  (  ■^  -  ). 

From  Figure  3,  we  see  that  tan  (  ^  )  =  r  A4>/Ax  .  Passing  to  the 

limit, 


"^x  4> 


(4) 


Substituting  the  stress -strain  relation,  ^  Tx  4*  ”  \  9  ’  obtain 

G  r-l^  =  T  rf,  •  (S) 

ax  x9 


Now  tlie  torque  exerted  on  the  cross  section  at  x  by  the  adjacent 
section  is  T  (x,  t ).  This  torque  is  a  result  of  the  shearing  stress,  ^ 
the  section,  so  that 


T(x,  t)  - 


J  G 


a  4> 
a  X 


(x,  t). 


From  (6),  for  a  shaft  of  constant  cross  sectional  area,  we  see  that 


over 

(6) 


a  T 
a  X 


(X.  t  )  =  J  G 


a^  4> 
a  x^ 


(x,  t  ) 


(7) 


Substituting  (7)  into  (3),  we  obtain 


5. 


(x,  t )  - 


2 


d  X 


(x.  t) 


L(x,  t) 
P  J 


(8) 


.  2  G 

where  c  =  —  • 

P 

Eq.  (8)  is  the  differential  equationof  motion  of  an  elastic  circular 
shaft  to  which  is  applied  a  time  dependent  torque,  L  (x,  t  )  ,  per  unit  length. 
It  is  valid  for  deformations  which  are  small  enough  to  insure  conformity  with 
the  behavior  assumed. 


6. 


Chapter  II 

FORMULATION  OF  THE  PROBLEM  AND  FORMAL  SOLUTION 


2. 1  We  wish  to  determine  the  torsional  motion  of  an  artillery  barrel  subjected 
to  the  time  dependent  torque  of  an  accelerating  projectile. 

Let  the  barrel  itself  be  represented  by  the  shaft  considered  in  Chapter  1. 
Its  equation  of  motion  then  is 


9^4) 


,2  9^<t> 


L  (x,  t ) 
P  J 


(8) 


9  t  9  X 

If  the  moment  imparted  to  the  barrel  by  the  projectile  is  considered  to  be  a  line 
moment  and  if  the  distance  of  the  projectile  from  the  left  end  is  designated  as 
^  ^  (t)  ,  then  the  following  equation  may  be  used  to  describe  the  physical 

system; 


2  2 
9  4)  2  9  4> 

- rr-  -  C  - IT- 


-ff  -  6  (X,  e  ) 


9  t  9  X 

where  6  (x,  ^  )  is  the  Dirac  Delta  function  having  the  property: 


(9) 


j 

0  for 

X 

Xf  = 

^2  - 

e 

e )  dx  == 

< 

1  for 

VII 

X 

VII 

4-01 

^2 

0  for 

e  s 

X 

Xf  = 

^2 

> 

and  M  (t)  is  the  time  dependent  moment  upon  the  shaft. 

We  now  obtain  an  expression  for  the  moment  M  (t):  The  projectile 
(assumed  rigid)  at  any  time  t  has  rotated  through  a  total  angle  0  0  (5,  t  ). 

This  angle  is  composed  of  the  angle  4>  due  to  the  deformation  of  the  barrel  and 
the  angle  a  due  to  the  helical  rifling  of  the  barrel.  Hence,  the  total  angle  through 
which  the  projectile  has  rotated  is 


7. 


e  =  4>(e.  t)  +  ke (t) 


(10) 


where  k  is  a  constant  dependent  upon  the  helix  angle,  so  that 


M(t)  =  1  6  =  I  [  ke  +  *<^(e,  t)] 


(11) 


where  1  is  the  mass  moment  of  inertia  of  the  projectile  and  the  dots  denote 
differentiation  with  respect  to  time. 

Eq.  (9)  may  now  be  written 

^  [  1^  V  +  4^  (S,  t )  ]  6  (x,  e ) .  (12) 

2.  2  Formulation  for  Barrel  Fixed  at  One  End  and  Otherwise  Free. 

Consider  Eq.  (8)  with  L  (x,  t  )  =  0  as  applied  to  Figure  4: 


=  0 . 
a  t 


a  X 


(13) 


% 


8. 


The  totality  of  motions  of  a  continuous  vibrating  system  is  known  when 
the  eigenvibrations  are  known.  1  j 

We  consider  free  vibrations  for  which  the  angular  displacement,  4)  , 

can  be  expressed  as  the  product  of  a  factor  X  (x),  depending  only  upon  position, 
and  a  factor  q  (t),  depending  only  upon  time. 

Substituting  4>  =  X  (x)  q  (t)  into  (13),  we  obtain 


1  q  (t)  ^  X"  (x) 

q  (t)  X  (x) 


(14) 


where  the  primes  denote  differentiation  with  respect  to  x  . 

The  left  side  of  (14)  depends  only  upon  time,  whereas  the  right  side 

depends  only  upon  position.  Therefore,  each  must  be  equal  to  the  same  constant 
2 

-X  ,  and  we  obtain  from  (14) 

X"  (x)  +  X  (x)  =  0 

..2  2  (15) 

q  (t)  +  X  c  q  (t)  =  0  . 


Solving  the  first  of  (15)  , 

X  (x)  =  A  sin  X  X  +  B  cos  X  x  .  (16) 

Now  the  boundary  conditions  are 


(a)  4)  (0,  t  )  =  0 

<b)  II-  <  £  .  t)=0. 


The  numbers  in  brackets,  unless  otherwise  designated,  refer  to  the  bibliography. 
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Condition  (a)  states  that  there  is  no  displacement  at  the  fixed  end.  Condition  (b) 
states  that  the  moment  at  the  free  end  is  zero. 

From  (a)  we  obtain 


B  =0. 


From  (b)  we  see  that  A  X  cos  (  X  1)  =  0.  But  if  A  =  0,  the  solution 

would  be  trivial.  Therefore  X  £  =  (2n  +  l)  .|  (n  =  0,  1,  2, . , ,  )  must  be 

taken  as  the  solution  in  order  that  the  boundary  conditions  be  fulfilled.  Hence 

the  solution  to  (13)  subject  to  conditions  (a)  and  (b)  above  is 

oo 


4>(x,  t)  = 


sin 


X 

n 


X  <t) 


(17) 


n=0 


Where  . 

The  functions  sin  X^^  x  are  the  eigenfunctions  of  the  differential  equation 
(13)  under  the  boundary  conditions  X  (0)  =  0  and  X'  (  £.  )  =  0. 


Orthogonality  of  Eigenfunctions  (x) 


We  shall  show  that  the  set  of  eigenfunctions  X^  (x)  is  an  orthogonal  set  over 

[o,  £]  ,  i.e.  , 

X  dx  =  0  for  m  ^  n 
m  ' 

To  this  end,  we  multiply  the  first  of  Eqs.  (15)  by  X^  (x)  ,  then  by  X^  (x)  and 
subtract  the  two  results.  This  yields 


X  X" 
n  m 


-X  X"  +(X' 
m  n 


m 


)  X. 


X  =0 
m 


(18) 


Integrating  (18)  over  the  length  f  , 


10. 


(X^  -  X^  )  XX  dx  +  (X  X'  -  X^  X'  )  dx  =  0  .  (19) 

'm  n'/nm  /dx'n  m  m  n 

Since  X  (0)  =  X  (0)  =  0,  X’^  (  >  =  X*  (  jt  )  =  0,  and  X^^  ^  X^  ,  (19)  states 


X  X  dx  =  0  for  n  ^  m  . 

I  m  n 

Jo 

Hence  the  X^^  (x)  =  sin  X^  x  subject  to  the  conditions  X^  (0)  =  0  and  X’^^  (  £  )  =  0 
form  an  orthogonal  set  of  functions  over  the  interval  ^0,  .  This  set  of 

functions  may  be  normalized  by  requiring  that 

[a^  sin^  X  X  dx  =  1  .  (20) 

/  n  n 

■^0  _  . 

As  a  consequence  of  (20),  ~”\/j  * 

The  trigonometric  functions  sin  X^  x  form  a  complete  set  on  the 
interval  j^O,  ,  and  therefore  any  piecewise  smooth  function  on  that 

interval  may  be  expanded  in  terms  of  sin  X^  x  such  that  the  series  expansion 
converges  uniformly  to  the  function.  (See  for  example,  [  ■*■  ] 

We  now  show  that  we  can  express  the  Dirac  Delta  function  in  terms  of 
the  eigenfunctions.  Let  y  =  p  (x)  be  a  continuous  function  such  that  p  (x)  =0 
unless  ^-e  =  x^  and  p  (x)  =  constant  =  for  x  in  that  interval. 

This  situation  is  depicted  in  Figure  5. 


Figure  5 
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We  can  expand  p  (x)  in  terms  of  sin  x  : 


P  (x) 


oo 

=1. 


sin  X  X 
n 


2  I  n=0  2  sin  X  ^  sin  X  e 

where  /  P  (x)  sin  X^^  x  dx  =  - ^ - - — 


0 


OO 


Therefore  p  (x)  =  j 

A 

<1 


»  X  sinX  ^  sin  X  e  sin  X  x 
^  \  n  n  n 


X  e 
n 


and  ior\[^  sin  X^  x  =  (x) 


n  n 


oo 


lim  /  , 

—  0 


lim 


n=0 


(i)  X  (x)  -  . 

n  '  '  X_  e 


sin  X  e 

"  =  6  (X.  e ) 


oo 


Xn  (e)  (X) 


(21) 


Utilizing  the  expressions  (17)  and  (21)  in  Eq.  (9),  we  obtain 
oo  oo 

V[x_,(x)qn-c2x-„(x)qJ 


n=0  2 

But  X"  (x)  =  -X  X  (x)  ,  so  that  (22)  becomes 
n  '  n  n  '  ' 


(X)  X^  (g) 


(22) 


niO 


oo 

I 

n=0  " 


2  M  (t)  X  (e) 

q  +  CO  q  -  - - 

^n  n  p  J 


X  (x) 
n  '  ' 


=  0  . 


2  2  2 
where  co  =  c  X 

n  n 

Because  the  functions  X^  (x)  form  a  complete  set  on  the  interval 

each  of  the  coefficients  in  (23)  must  be  zero.  Hence, 

2  M(t)X„(e) 


q  (t)  +  CO  „  q  (t)  = 
m  n  ^n 


P  J 


(23) 


[o.^] 


(24) 


qn  (t). 


must  be  solved  for  the 


12. 


2.  3  Formulation  for  Barrel  Free  at  Both  Ends  With  Two  Elastic  Torsion  Supports 
Along  its  Length. 

This  system  is  depicted  in  F  igure  6. 


The  differential  equation  for  this  case  is  the  same  as  that  previously  considered, 
namely 


,2  <!> 


a  X 


M(t) 
P  J 


6  (X,  e )  . 


The  supports  at  a  and  b  resist  the  angular  displacement  of  the  barrel  such  that 


where 


across 


T^  =  4>  (a,  t  ) 

and  02  are  constants. 

The  boundary  conditions  which  must  be  satisfied  are; 
the  supports  must  be  continuous.  This  condition  yields 


(25) 


(1)  the  moments 


J  G  X’^  (a‘)  -  J  G  X'^  (a+  )  =  -a^  X^  (a) 

J  G  X'^  (b‘  )  -  J  G  X'^  (b+)  =  -a^  X^  (b)  ; 


(26) 


13. 


(2)  4>  must  be  continuous  at  a  and  b  .  Therefore 


Xn(a-)=Xn(a^> 


(b' )  =  (b^  )  ; 


(27) 


(3)  the  ends  must  be  without  moment  so  that, 


X’n(O)  =  X'^(  1  )  =  0 


(28) 


We  know  that  X^^  (x)  must  satisfy,  in  addition  to  the  six  boundary 
conditions  given  by  (26),  (27),  and  (28), 

<*>  +  \  -  0 

In  each  of  the  three  regions.  The  differential  equation  can  be  satisfied  by  the 
functions  X^^  (x),  defined,  for  convenience,  by 

X  (x)  =  u  (x)  =  A  cos\  (a-x)  +  B„  sin  (a-x)  for  OS  x^  a  (29) 

n  n  n  n  n  n 

Xn  (x)  =  Ti^  (x)  =  cos  X  +  sin  for  a  ^  x  =  b  (30) 

X  (x)  =  t  (x)  =  E  cos  \  (x-b)  +  F  sin  \  (x-b)  for  b=x=  £.(31) 

n''^n  n  n  n  n 

By  (26)  and  (27)  we  find  that 


JG\  B  +JGX  (D  cos  X  a-C  sinX  a)=A  a, 
■’  nn-*  n  n  n  n  n  n  l 


=  a,  (C  cos  X  a  +  D  sin  X„  a  ) 
1  '  n  n  n  n 


(32) 


and 
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-J  G  F  +  J  G  X  (D  cos  X  b  -  C  sin  X  b  )  =  -E  a„ 
•’  nn-’  n'n  n  n  n  n  2 


=  -a,y  (C  cos  X  b  +  D  sin  X  b  )  , 
z  n  "  ^  " 


n  n 


so  that 


A  =  C  cos  X  a  +  D  sin  X  a 
n  n  n  n  n 

=  C  cos  X„  b  -f  D  sin  X  b 
n  n  n  n  n 

A  a. 

B  =  ttS-t —  +  C  sin  X  a  -  D  cos  X  a 
nJGX^  n  n  n  n 

+E 

F  =  ,  . - -  C  sin  X  b  +  D  cos  X  b  . 

nJGXj^  n  n  n  n 

By  condition  (28)  we  have 

A  sin  (X  a  )  -  B  cos  (X  a  )  =  0 
n  n  .  n  n 

-E  sin  X  (i  -b  )  +  F  cos  X  (  f  -  b)  =  0  . 
n  n  n  n 


(33) 

(34) 

(35) 

(36) 

(37) 

(38) 

(39) 


Substituting  expressions  (34),  (35),  (36)  and  (37)  into  (38)  and  (39)  we 
obtain  two  simultaneous,  homogeneous,  algebraic  equations  in  and  . 

In  order  that  the  equations  have  a  non-trivial  solution,  the  determinant  of  the 
coefficients  must  be  equal  to  zero.  This  condition  produces  the  following 
equation  in  X^  ; 

2  2 

(J  G)  X  sin  X^  j£  +  J  G  X^  ^  cos  X^  a  cos  X^  (  £  -  a  ) 

+  o,  cos  X  b  cos  X„  (  f  -  b) I  +  a,  a,  cos  X  a  cosX  ( £  -b)  sin  X  (b-a) 
z  n  n'  'J  iz  n  n  '  n  ' 

=  0 


(40) 
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The  roots  of  (40)  are  the  eigenvalues  for  the  problem  of  the  barrel  having  two 
elastic  torsion  supports  along  its  length. 

The  eigenvalues  having  been  obtained  by  Eq.  (40),  we  solve  the 
homogeneous  equations  in  the  following  form: 


C  =  c  4> 
n  n  ^n 


D„  =  d  ill 
n  n  ^n 


We  are  at  liberty  to  choose  one  of  the  c^  or  d^  equal  to  unity  and  then 

determine  the  other  by  either  of  Eqs.  (38)  or  (39).  The  tp  however,  sire 

n 

completely  arbitrary.  We  may,  therefore,  choose  them  so  that  the  functions 
(x)  are  normalized,  i.  e. ,  so  that 


I  X  j^(x)  dx  =  1  . 

®  r  1 

By  the  method  of  Rayleigh  2  ,  Eq.  (41)  may  be  written 

-£  r  ^  ^  £. 

x(-^)^  +  xX^  - -i-  X  X’ 

n 


X\(x)dx=  i 


'0 


=  1 


0 


(41) 


(42) 


Therefore,  for  the  problem  at  hand,  Eq.  (42)  demands  that  the  ip^  be  chosen  such 


that 


+  b 


<b)  -  (b) 


+  £  £"„(£)+ 


(43) 


We  seek  to  show  that  the  set  of  functions  X^  (x)  defined  by  Eqs.  (29), 
(30),  and  (31)  is  orthogonal.  To  this  end  we  write  the  differential  equation  for 
two  individual  terms  of  the  set: 
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X’’  (X)  +  X  „  X  (X)  =  0 
n  n  n 


X"  (x)+x"  X_(x)  =0 
m  mm 


(44) 


If  we  multiply  the  first  of  (44)  by  X^  (x)  and  the  second  by  X^^  (x),  subtract  the 
two,  and  integrate  the  difference  from  0  to  jg  ,  we  obtain 


2  .2 
m  n 


/  ^rry  (*)  ^r.  (*)  ^  + 

X  X' 

-  X  X’ 

I  m  '  ^  n  '  ' 

n  m 

m  n 

0 

L  J 

=  0  (45) 


Expanding  the  second  term  of  (45),  we  obtain  by  definition  (29),  (30),  and  (31), 


X„  X'  -  X^  X’ 
n  m  m  n 


=  (a)  (a)  -  AX^(a)  (a)  -  (0)  (0) 


By  (28), 


(46) 


■^n  ^'m  +^m 

+  fn  S  m<^>  -  fm<^)  £  n<^)  -  fn<'^)  +^m  ^’n 


M'n(0)=  for  all  n 


(47) 


By  (27), 


^n  for  all  n. 


(48) 


Hence 


XX'  -  X  X' 
n  m  m  n 


=  (a) 


M'  (a)  -  Ti’  (a) 


-li  (a) 


-  h  (b) 
m 


1*’  (a)  -  Ti'  (a) 
n  n 


m 


+  hn  (b) 


m 


V^(b)  -  f^b) 


Vn(b)  -f„(b) 


(49) 


% 
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By  condition  (26) 


(a) 

’’ %  W  forall  m.  (5 

Therefore  the  aecond  term  In  (45)  is  zero,  and,  since  X„  in  general  Is  different 


from  for  m  n  , 


/  (x)  d  X  -  0  for  m  n  . 

-^0 

Hence  the  (x)  defined  by  (29)  through  (31)  subject  to  (26)  through  (28)  form 
an  orthogonal  set  of  functions  on  the  interval  ^0,  fj 

The  formal  solution  to  Eq.  (24).  to  which  both  systems  previously 

discussed  reduce,  is  given  by  Duhamel's  integral  , 

_f 


qn(^)  =  -i — 
a  P  J  w 


I  +■-!>'[  e  (f).  t']  I  sin  u,^  (t-f)  d  f 


a  sin  w  t  +  b  cos  co  t 
n  n  n  n 


If  the  projectile  begins  from  rest  at  the  left  end  of  the  barrel,  then 
\  \  °  <5i)-  If  we  multiply  (51)  by  (^)  and  sum  over  n  ,  we 

obtain  the  expression  for  the  deflection  at  the  line  moment: 

oo 

^n  «>%<'>' p-J  /  ke(f)  +  +  [e<0. 

n=0  "^0  '  /  '  “n  / 


=  «<>(€.  t)  =  <(.^(0 


We  may  reduce  the  integrodifferential  equation  (52)  to  an  integral  equation  by 


noting  that 


t 


18. 


<)>o(t)  =  /  (t  -  f)  d-*. 


Relation  (53)  may  be  verified  by  direct  integpration: 

/  (t  -  f)  4;^  (f)  df  =  -t  4.Q  (0)  +  <|>^  (t)  -  45  (0)  . 

By  the  condition  that  4>  (x.  0)  =  4.  (x,  0)  =0,  (54)  reduces  to  4>q  (t). 
Substituting  (53)  into  (52),  we  obtain 

^[ke  (t')  K(t,r)  -  (t-f)]  dt*=  J[(t-t')  -  K(t,  t-)]  ■<ij(.t')dt' 

r  1  r  1 


where 


M.,.) 

pJA  « 


sin  (t  -  t') 


(53) 


(54) 


(55) 


The  solution  of  (55)  presents  foimidable  mathematical  difficulties.  An 
approximate  technique  will  therefore  be  employed  in  order  to  gain  some  knowledge 
of  the  barrel  behavior. 


Chapter  III 


19. 


APPROXIMATE  SOLUTION  AND  NUMERICAL  RESULTS 

By  taking  only  the  first  few  terms  in  the  series  given  by  Eq.  (17),  the 
solution  of  (24)  can  be  reduced  to  that  of  a  system  of  simultaneous  ordinary 
differential  equations.  With  this  end  in  view,  let 


«*>  (X.  t )  =  Xq  (X)  q^  (t)  +  X^  (X)  q^  (t) 

Then  Eqs.  (24)  become: 

%<=•  -  Qi/Jo/!,  - 

Pr\  fin  ^ 


o^o^o  +  ^  ^ 


and 


(56) 


(57) 


q'l  (c*  q-Q  /3q  -  2  q^  -  2  q^ 

-  q^  p^  p^^  +  c*  %Pi  Pq  =  k  *e  . 

Multiply  (57)  by  (c*  -  p^^  )  and  (58)  by  p^  p^  and  add  to  obtain 

•  '  P  0  -  =  2  Qq  /^o  ^0  ^  ^1  ^0  ^1 

-(c*  CO  Q  -  /3q  p'^ ^  )  q^  -  (p^  ^  -p^  )q^+  kp^ 

Similarly,  multiply  (57)  by  p^  p^  and  (58)  by  (c*  -  p^^  )  and  add  the 
two  to  get 


(58) 


(59) 


20. 


Ql  (c*  -  -P^^)=2q^  +2^p^p^ 


-  (c*  -  Pj^  p\  -  p\  )  q^  -  (p^  p^  <^^0  -  Pl'p'o^  %  I'. 


(60) 


where 


=  £J. 

I 


and 


P  =  X  (g) 


Eqs.  (51)  and  (52)  can  therefore  be  solved  for  the  two  generalized  coordinates, 
Qq  (t)  and  q^  (t)  . 

The  solutioiBto  these  equations  were  obtained  in  graphic  form  by  means 

of  an  analog  computer.  By  comparing  orders  of  magnitude  of  the  various 

terms  in  equations  (59)  and  (60)  and  due  to  the  limitations  of  the  computer, 

all  terms  involving  derivatives  of  the  P’s  were  neglected  in  these  equations. 

2  2 

Also  the  term  (P  q  +  p  ^  )  was  neglected  as  compared  with  c*.  With  these 


simplifications,  Eqs.  (59)  and  (60)  become  respectively 


V  _  “^0  ^^1  2  „  ^0  ^1  "^1 

9o  p  90  "  0  ^0  c* 


91+ 


"  ^0  V 


(61) 


9,  = 


2  fl2 
‘^1^0 


^1  1^1 


^0  ^1  “  0 


^p,  .. 

9o+  -p-  ?  . 


(62) 


It  should  be  noted  that  the  above  simplifications  have  been  made  primarily 
due  to  computer  limitations;  some  of  the  terms  neglected  are  comparable  in 
magnitude  to  some  retained. 


Numerical  Results 

Some  typical  artillery  data  used  in  determining  the  coefficients  in 
Eqs.  (61)  and  (62)  are  given  in  Table  I.  The  calculated  values  of  the  problem 
parameters  are  tabulated  in  Table  11. 


TABLE  I 


Typical  Artillery  Data 


Weight  of  gun  tube 

Length  of  tube 

Projectile  travel 

Weight  of  standard  H.  E.  round 

Inner  tube  diameter 

Le  Due's  Velocity  Formula 

Maximum  muzzle  velocity  for  zone  7 

Helix  angle 


1058  lb 
7.  75  ft 
6.  8  ft 
33  lbs 


4.  13  inches 


V  = 


1784  g 
1.0277 


1550  ft/ sec 


8°  56’ 


TABLE  II 

Problem  Parameters 


Xq  =  tt/2  £  =0.231  ft'-^ 
Xl  =  =  0.693  ft'^ 

G  (steel  )  =  12  X  10^  psi 
P  (steel  )  =  488.  8  Ibm/ft^ 
^0  ~  2465  rad/sec 


cjj=\y  —  Xj^  =  7394  rad/sec 
,  _  2  tan  8°  56'  j/,.. 

*=  -  -ni  t,- 


sin  Xq  I  =  .54  sin  X^  g  ft 
sin  X^  I  =  .  54  sin  X^  g  ft 


j. 

2 

_1 

2 


c*  (for  shell  weight  of  33  lbs.)  =  20.  67  ft 


22. 


In  calculating  the  value  of  the  polar  moment  of  inertia  of  the  tube  J, 
the  tube  weight  and  the  length  of  7.  75  ft  were  used;  in  all  other  calculations, 
however,  the  tube  length  and  the  projectile  travel  were  taken  to  be  the  same 
value  of  6.  8  ft. 

Figures  7  through  10  show  the  variations  of  and  with  time  and 
distance  of  the  projectile  down  the  barrel  for  several  different  values  of  the 
parameter  c*.  In  Figure  11  is  given  the  variations  of  g,  g  and  g  with  time. 
Figures  12  and  13  show  the  variation  of  the  q’s  and  q’s  with  c*  at  the  instant 
the  projectile  leaves  the  barrel.  We  note  that  the  angular  displacement  and  the 
angular  velocity  of  any  section  of  the  barrel  are  relatively  small  at  this  time. 

In  Figure  14  is  plotted  the  angular  displacement  (j>  (x,  t )  as  a  function  of  x  for 
the  nominal  value  of  c*  equal  to  20.  67  and  for  three  values  of  time;  tg  ,  the  time 
at  which  the  projectile  leaves  the  barrel  ;  0.  6  tQ  ,  the  approximate  time  at 
which  the  sum-  q^  +  q^  reaches  a  maximum;  and  0.  92  t^,  the  time  at  which  q^ 
reaches  a  minimum. 


# 


Figure 


Figure 


F igure  10 


8.7 


Figu 
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APPENDIX  A 


In  this  appendix  is  presented  the  analog  computer  program  used  in 

solving  Eqs.  (61)  and  (62)  ,  which  are  rewritten  here  for  convenience. 

2 


^0  = 


qi  = 


0^  1 


2  «2 
‘"1^0 


Qi  -  1  - 


^0  “  1 


%  -  0 


^1  + 


^0  + 


(61) 


(62) 


To  increase  the  computer  time  to  approximately  40  sec. ,  a  time  scale  factor 
of  4000  was  chosen.  For  the  nominal  value  of  c*  equal  to  20.  67  ,  a 
preliminary  calculation  determined  that  it  would  be  convenient  to  scale 
q^  and  q^^  by  factors  of  10^  and  5  x  10^  respectively.  The  voltage  and 
time  scaled  equations  then  are 

6  “^0  5 

)  qo  =  ( —A  ^  )  lo^qn  -*  ( 


and 


10 


a  c” 


0  V  ,^,5 

10  qo 


,  2  Pq  5  •  iO^  k  /3 

— -ac*  -  >  5x10^  q^+( - 10^ 


6  /  ^  ^  I  ^  0  5.  2w 

<  — ac*  )  5  X  10^  q^  -  (  -^)  5  x  10^  q^ 


-( 


a  c’ 


•)  10  qp  +  ( 


)  10"^', 


(63) 


(64) 


where  a  =  16  x  10  and  the  primes  denote  differentiation  with  respect 
to  T  =  4000  t  . 

The  function  ^  was  chosen  to  approximate  an  experimental  pressure 
curve.  The  analytical  expression  chosen  was 


33. 

"^  =  A  e  ( 1  -cos p t  )  .  ^ 

where  p  =  2ttx10  and  0=  t  ^.01.  The  constant  A  was  determined  by 
requiring  that  the  maximum  muzzle  velocity  be  1550  ft/sec  and  a  was  found 
by  requiring  that  the  peak  acceleration  occur  at  .  38  p  T*  where  T*  is  the 
period.  Utilizing  these  conditions 


and 


A  =  4.  83  X  10^  ft/sec^ 
a  =  .393  . 


Eq.  (65)  was  put  into  suitable  form  for  programming  by  integrating 
twice  and  eliminating  the  trigonometric  terms.  The  two  constants  of  integration 
were  determined  from  the  conditions  that  e  (0)  =  g  (0)  =  0.  An  appropriate 
voltage  scale  factor  for  this  equation  was  found  to  be  10^  .  The  equation 
which  was  put  on  the  computer  is  as  follows: 


10  e"  = 


-  ai 


e* 


2  2 
P  (1  +  a  )  e 

16  X  10^ 


- 1 — T 

16  X  10  a 


(e 


-apT 

4000 


-  1) 


+ 


Apr 

64  X  10^  a 


(66) 


Substituting  the  numerical  values  of  a,  p,  and  A  into  Eq.  (66),  we  obtain 

10^  e"  =  -  123  e'  -  28.5e  +  198.3(  e  --0618  T  .  i  +  ^2  t  . 

The  computer  block  diagram  for  the  solution  of  Eqs.  (63)  and  (64)  is 
shown  in  Figure  (A.  1), 
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